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■ ABSTRACT 
D ■ 

. Context. We investigate the growth of hydromagnetic waves driven by streaming cosmic rays in the precursor environment of a su- 

' pernova remnant shock. 

' Aims. It is known that transverse waves propagating parallel to the mean magnetic field are unstable to anisotropies in the cosmic 

ray distribution, and may provide a mechanism to substantially amplify the ambient magnetic field. We quantify the extent to which 

. temperature and ionisation fractions modify this picture. 

Methods. Using a kinetic description of the plasma we derive the dispersion relation for a collisionless thermal plasma with a stream- 
ing cosmic ray current. Fluid equations are then used to discuss the effects of neutral-ion collisions. 

Results. We calculate the extent to which the environment into which the cosmic rays propagate influences the growth of the magnetic 
^ ■ field, and determines the range of possible growth rates. 

5— ( ' Conclusions. If the cosmic ray acceleration is efficient, we find that very large neutral fractions are required to stabilise the growth 

' of the non-resonant mode. For typical supernova parameters in our galaxy, thermal effects do not significantly alter the growth rates. 

, For weakly driven modes, ion-neutral damping can dominate over the instability at more modest ionisation fractions. In the case of a 

1 ' . supernova shock interacting with a molecular clouds, such as in RX J1713.7-3946, with high density and low ionisation, the modes 
' can be rapidly damped. 

^ ' Key words, plasmas - instabilities - cosmic rays - supernova remnants - magnetic fields 

cn : 
^ . 

1. Introduction 

rn 



There is clear evidence from radio, X-ray and gamma ray observations that electrons, and probably protons, are accelerated in 
f-*) supernova remnants (SNR). First order Fermi acceleration is believed to be the most likely mechanism for this acceleration, as it 
f — provides a natural explanation for the spectral shape (for a recent review see Hillas 2005). Cosmic rays up to the ankle are believed 
to have their origin in our galaxy. However, the standard picture of Fermi acc eleration must be stretched to its limit in order to 
;V accelerate protons to even the knee of the spectrum at ~ 3 x lO'^ eV (Lagag e"& Cesarskvl [l983.) . Various methods have been 
. !^ suggested to solve t his problem in the qua silinear framework, but amplification of the ambient magnetic field probably provides the 
best possibility (e.g. lKirk & Dendvll200ll) . The presence of bright X-ray rims near the shocks of many young supernova remnants 
has been interpreted as observational evidence of field amplification ( Vink & Laming 2003; Volk et al. 2005; Ballet 2006), but more 
detailed comparisons of these data with radio observations are necessary to confirm this result (Katz & Waxman, 2007). 

Recently an attempt to derive a non -linear theor y of magnetic field amplification was made by Bell & Lucek (.200 L) . and sup- 
ported by hybrid numerical simulations dLucek & Be ll. 2000). It was shown that the cosmic ray pressure could wind up t he frozen-in 
magnetic field to values in excess of the typically assumed self-quenching level of SB/Bq < 1. Building on these papers. Bell' ( 2004|) 
reanalysed the work of Ach terberg ( 1983), and found a strong non-resonant mode that had been previously overlooked. Numerical 
MHD simulations inlil 12004, 2005) verified that the streaming of cosmic rays can amplify the magnetic field at short length- 
scales to values greater than the ambient field. The instability is driven at short wavelengths due to the thermal plasma attempting 
to neutralise the effects of the streaming cosmic rays. This induces a current in the thermal plasma. The cosmic rays essentially 
behave as unmagnetised particles and do not move cohe rently with magnetised thermal background particles. The same instability 
has also been seen, using particle in cell simulations, (.Niemiec & Pohl. .2007) . These authors observe the non resonant instability 
with a considerably slower growth rate, but this may be due to the artificially low proton to electron mass ratio {nip/in^. -10) used 
in their simulations. This considerably reduces the range over which modes can be unstable. More realistic mass ratios are required 
to understand the true nonlinear development of this instability. 

The upstream medium through which the shock-accelerated particles stream, w ill, in general, be in completely ionised. Limits 
on the acceleration of cosmic rays in these environments have been calculated bv iDrurv et al.l (Il996l) . including the damping of 
resonantly excited Alfven waves due to friction with the neutral particles. This will also play a significant role in molecular clouds 
where the ionisation fractions are even lower and densities much higher The propagation of Alfven waves in molecular clouds has 
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been studied by iTagger etaTI (Il995h . Extending the work of iKulsrud & Pearcd (Il969h and lZweibel & Shulll (Il982h . they demon- 
strated that the ponderomotive force of the waves can have a significant effect on the structure and resulting magnetic field within 

such clouds. 

Recently. iBvkov & ToptvginI (|2005|) . have investigated a similar instability to that in Bell (2004), driven by a neutralising return 
current in the thermal background plasma, in a partially ionised medium. A single fluid MHD description of the background particles 
with a generalised Ohm's law is used, and the cosmic ray current is found by solving the kinetic transport equation. They conclude 
that a non-resonant mode with a rapid growth rate remains. However, as we show in section 3, for high frequency waves the friction 
between the charged species and the neutrals is not sufficient to couple their motions, and a single fluid description of the plasma is 
no longer appropriate. 

In the present work, we investiga t e vario us conditions that may reduce the growth of the non-resonant instability. In section 2 
we follow the analysis of lAchterberd (Il983h . neglecting collisions, and quantify the extent to which thermal effects can reduce the 
growth of the instability. In section 3 we introduce a two fluid description of the plasma which includes the collisions between the 
different species. The neutrals are coupled to the ions through the ion-neutral friction, and the cosmic rays are coupled to the MHD 
fluid through Maxwell's equations. Section 4 applies our results to SN1006 and RX J1713. 7-3946. We conclude with a discussion 
about some of the consequences of the instability in various conditions. 



2. The kinetic non-resonant current driven instability 

The linear dispersion relation for ckcularly polarised transverse waves propagating parallel to the zeroth order magnetic field is 
— -l = yxAk,(o), (1) 



where the summation is over each species with corresponding charge qs, cyclotron frequency Wcs = ^sBo/nisC and plasma frequency 
= Ane^njnis. The susceptibility 5 for each component of a weakly ionised plasma is determined by integrating the Boltzmann 

transport equation along the unperturbed trajectories about the zeroth order field (e.g. Krall & Trivelpiece 1973). 

Previ ous work on the eff'ect of streaming cosmic rays using a kine tic approach has assumed the background plasma to be colli- 

sionless (lAchterbergi[T983l:lReyille etaill2006tlBlasi & AmatoLl2007h . For a three-species plasma with proton electron background 

and proton cosmic ray component, Achterberg determined the background susceptibility for low frequency waves, tbs <s: a>d < \(i>ce\, 

to be 



.2 k^Vl_ 
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(2) 



where Cjs - uj- ku^ is the Doppler shifted frequency with the group drift speed of each species, \a - Bq/ V47!'ni'«i is the Alfven 
velocity and V^^ - ksTi/mi the ion thermal velocity. For ease of notation we have also introduced the parameter e to describe the 
polarisation with e = +1 for right-handed waves and e = -1 for left-handed waves (oj > 0). The requirements of zero net current 
and quasineutrality in the shock precursor, where the proton cosmic ray component is anisotropic, induces a compensating drift in 
the background plasma 7^'^ = n„ (ua - oj/k), where n„ and u^- are the cosmic ray number density and bulk speed, respectively. We 
note that n^i is a function of distance from the shock front, since only high energy cosmic rays are expected to penetrate far into the 
precursor 

Using the kinetic treatment of Bell for the cosmic ray current, the resulting dispersion relation takes the form 



lk^V}.\ 
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where, 

as appropriate in the rest frame of the upstream plasma in the precursor of a shock front, we have identified with the incoming 
plasma speed Vj. We define - p^inc/eBo the gyroradius of the minimum energy cosmic ray and consider only waves with ^ > 0. 
The function cr represents the cosmic ray susceptibility normalised to the return current J^^ as described in Eq (12) of lBelll ( |2004|) . 
We extend this result for arbitrary power law distributions / oc p""*^*, between p^am and pm-^x- For q < I 
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where 2F\{a, b, c; z) is the Gauss hypergeometric function. 

In (O we have also introduced a dimensionless parameter characterising the strength of the driving term 



(4) 
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In a shock precursor, we expect /j^in to increase with distance ahead of the shock front (e. g.. lEichleA 1 1 979t iBlasil l2002h . since 
higher energy particles have a larger mean free path. For ^ = 0, and neglecting a logarithmic term, this leads to «cr °c l/Pmin- Thus, 
^ is approximately constant, reflecting the fact that the energy density of the cosmic rays is the same in each decade of momentum. 
It is roughly given by Vj/c times the ratio of the cosmic ray energy density at the shock to the incoming flux of momentum riiniiV^. 
If SNR shocks are efficient in accelerating cosmic rays, one expects this ratio to be at least several per cent, (e.g., Volk et al., 19841). 
Hence, ^ ~ O.OlVs/c. For a more detailed discussion of this quantity see section 7 of Bell (2004). 

The dispersion relation ([3]) clearly permits waves travelling in directions parallel and antiparallel to the mean magnetic field for 
both left and right polarised waves. For the case of parallel streaming of the cosmic rays along the magnetic field, the waves that 
propagate parallel the field (w > 0) are unstable, and waves propagating antiparallel (w < 0) are damped. This can be understood 
physically, as the Lorentz force either acting centripetally or centrifugally on the field line, depending on the polarisation of the 
wave. 

We see from (|4|i that crikr^^) is a monotonically decreasing function of wavenumber for kr^^ > 1. The non-resonant mode 
emerges when cr is negligible, and the driving term dominates over the "Alfven term" that contains Va, i.e., when 1 <K kr^^ <sc 

^v2/v2. 

In the frame co-moving with the drifting thermal protons, the maximum growth rate of the non-resonant mode is 

Im(w) = * Wd, (6) 

2 Va 2 Va m 

which is independent of magnetic field strength. 

1 /2 

For a fixed f , the growth rate scales with the density as , but the mode grows on shorter length scales as rii increases, 
rendering it liable to thermal damping. Defining the dimensionless temperature - kgT jmiC^, for large k, and again neglecting cr, 
it follows from (l3]i, in the limit » m jc^, that the maximum growth rate is 

Im(w) ~ (n„/«i)'^' (jS?/0)'^' «ei, (7) 

where jSj = Vs/c. This is the non-relativistic equivalent of the result found bv lReville et alj ( l2006h . At high temperatures the thermal 
particles behave as though they were unmagnetised. For parameters appropriate for SNR, this effect is shown in Fig[T] In order for 
the non-resonant mode to leave the regime of linear growth before being overtaken by the shock front, i.e., before being advected 
over a distance of roughly crgm/Vj, one requires Im(<y) > v^/rgmC. The necessary condition for thermal effects to reduce the growth 
rate below this value is 

> ^2^. (8) 
miVs^ 

For typical SNR parameters this condition will only be satisfied for a very weak driving term , ^ <sc 1. From Fig.[T|it is clear that 
damping effects do not suppress the instability in the interstellar medium of our galaxy for typical SNR shock speeds of around 
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5000 km/s. Thermal effects are likely to pl ay a more significan t role for relativistic shocks, and may even provide a saturation 



Km/s. inermal errects are liKely to pl ay a more sigmncan 
mechanism for the current driven instabihty (iReville et al.L l2006l) 



3. The collisional non-resonant instability and ion-neutral friction effects 

Unhke the analysis of the previous section, the ideal M HD equation s implicitly include the effects of collisions, although they 
do not appear in the final equations. Thus the analysis of iBelll (|2004|) included the most important effects of Coulomb collisions. 
However, the media into which supernova shocks propagate are generally not completely ionised, and may even include weakly 
ionised molecular clouds. Consequen tly, we now extend t he ideal MHD equations to include col lisions between charged a nd neutral 
particles, following the approach of Tagg er et~aLl (1 1995b . who, in turn, extended the work of Kulsrud & Pearc3 (Il969l) . Whereas 
previous work involving ion-neutral friction has concentrated on shear- Alfven waves, we consider the non-resonant modes driven 
by the streaming cosmic rays. These propagate along the mean field and are, therefore, circularly polarised. We continue to use a 
kinetic description of the cosmic rays, which remain collisionless. Neglecting viscosity and Ohmic friction, the equations of motion 
for a frictionally coupled system of an MHD fluid and neutrals are given by 

dui j X B 

— = Vin(Ui - Un), (9) 

at piC 
and 

— = -Vni(Un - Ui), (10) 

at 

dTagger et al.L Il995h . where j is the current carried by the plasma (excluding the cosmic rays), Uj and Un are the velocities of the 
ionised and neutral components, respectively, and Vin and Vni are the momentum exchange frequencies. In writing these equations 
we have neglected all terms of order me//«i,n- The cosmic ray current is contained in Ampere's law V x B = (4;7r/c)(j 
The momen tum exchange frequencies a re related by PnVm - Pi^in^ where p„ and pi are the densities of the neutral and ionised 
components. lKulsrud & Cesarskvl (Il97ll) give the following approximation for plasmas in the temperature range 10^ K< T < 10^ K 

Vin ^8.9xlO-\„(^^)''' s-'. (11) 
The resulting dispersion relation is easily found to be 

aj^(l + J^]^kW + e(—k(cr-l). (12) 

In the high frequency limit, |w| » Vin, this reduces to a dispersion relation similar to (O, although the term representing thermal 
effects is absent. In the opposite limit |w| v^i we find 

2 

(13) 
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+ e^—k(cr - 1) 



For high collision frequencies, the neutral and ionised components are tied together and the effect of ion-neutral collisions is simply 
to increase the effective mass of the ions. 

For strongly driven, non-resonant modes (fT2l l reduces to 

,2 



(1 +Z)a? +Z(—k 



= 0, (14) 



where Z = pjpn = >'ni/i'in- For all physically relevant parameters (i.e., ^, Z > 0), this cubic polynomial in a> has three purely 
imaginary roots, two of which are damped modes, and the third the non-resonant growing mode. Thus, ion-neutral collisions are 
unable to stabilise the strongly driven mode, although they can affect its growth rate. 

The general analytic expressions for the growing mode are cumbersome, but various limiting cases yield interesting results. In 
the limit of low ionisation, Z » 1 one recovers the growth rate given in For Z 1 the growth rate is given by 

(15) 

If v^^ <K 4(vlk/rgjn collisions are too slow to compete with the driving, and the growth rate is again the same as ^ with only a 
small reduction due to collisions. However, in addition to the low frequency limit in which the components are tied together ST3[ . 
a new regime arises. If the collisional drag on the ions is sufficient to affect the driving, but not to couple the components, i.e., 
^fn '^i^sk/fgrn s> v^j, then we obtain a similar expression to that of lBvkov & Toptvginl ( 20051) . 




V,^ OJci ricr 

7=«^^— — (16) 

l^in^ gm ^in 

However, this applies only when y <K Vi,,. 

In Fig.|2]we illustrate the influence of ion-neutral collisions for SNR parameters. Thermal damping is negligible, but a weak 
dependence on temperature enters via Eq ( fTTT i. Using ( fT2l i. we plot the maximum growth rate, as a function of ionisation fraction 
Xi = njn, taking the driving term ^ = O.OlxiVs/c since the cosmic rays couple directly only to the ionised component. For higher 
density plasmas, a larger fraction of neutral particles is necessary to reduce the growth rate below the threshold value of v^/rg^c. 
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Fig. 2. Maximum growth rate as a function of ionisation fraction for total density n = + «„ in a H-H^ gas: 0. 1 cm ^, 10"* K (solid), 
1.0 cm-\ 10^ K (dash), 10.0 cm-\ 10^ K (dash-dot). We take ( = Q.QlxiV Jc with shock speed = 5000 km s"', and r^^ = 10''*cm. 
CO is in units of v^/rgmC. Once the comax ^ 1 the non-resonant instability is not effective in amplifying the magnetic field. 



4. Application to SNR 

X-ray observations of several young supernovae display bright rims immediately downstream of the outer shock. This has been 
modelled as the sy nchrotron coo ling of relativisti c elec t rons, which requi res a magnetic field substantially higher than in the sur- 
rounding medium dVink & Lami ng. 2003; Ballel 120061; IVolk et all l2005l) . In the case of SN 1006 the magnetic field at the shock 
front is estimated to be of the order 100//G. One way in which field amplification may take place is if the non-resonant instability 
discussed above occurs upstream of the shock front. Since this instability depends on strong driving by cosmic rays, it cannot boot- 
strap from an initial state in which they are absent. Nevertheless, we can check its consistency, by requiring that, given a substantial 
cosmic ray pressure at the shock, the growth rate is sufficiently rapid for it to enter the nonlinear regime before being overtaken 
by the shock. Thus, the physical conditions far upstream (B, n, x-,, T) together with the shock speed v, and cosmic ray intensity 
UalP'i^s^ - 0.1, should combine to yield a growth rate in excess of Vs^/(rgmc). 

In the case of SN1006, the sh ock velocity is a pproximately Vj ~ 2900 km s"' the density is 0.05 < n < 0.25 cm"^, and the 
neutral fraction is small, ^ 0.9 jRavmond et al.L [2007). Taking «i = 0.25 cm"^ gives ^Vj^/va^ ~ 10. The low density and shock 
speed, imply that the non-resonant mode is not very strongly driven, but the driving is nevertheless strong compared to the collision 
frequency. Very close to the shock, the growth time of the field is on the order of years, while further from the shock the growth rate 
decreases quite dramatically, as the number density of cosmic rays decreases. For the parameters given above the maximum growth 
rate is ~ lOvs^/rgmC. The non-resonant mode may still amplify the field into the nonlinear regime. 

The low neutral fraction in SN1006 means that ion-neutral collisions are not likely to dramatically reduce the instability. The 
interaction of a supernova blast wave with a molecular cloud, such as that observed close to the north-western rim of RX J 17 13.7- 
3946 would have a much more significant effect on the growth. Although the exact parameters are uncertain, mo lecular clouds ar e 
generally clumpy with interclump densities in the range 5-25 cm" and ionisation fractions not larger than 10% (IChevalieii[T999l) . 
Adopting the values in .Malkov et al. (2005) T = 10^ K, n = 23 cm"^, Xi - 0.01, Vs - lO^cm s"' the collision frequency is larger 
than the growth rate. The non-resonant driving term is unable to dominate over the Alfven term and growth is driven resonantly by 
the cosmic rays themselves, at very long wavelengths kr„jn <K 1 . At shorter length scales, the waves are rapidly ion-neutral damped 
Alfven waves. The shortest growth timescale is on the order of ~ 10^ yrs. This value varies with choice of rgm, temperature and 
density, but for high density and low ionisation fractions, the growth timescale is longer than the free expansion phase of a supernova 
remnant, or the lifetime of the wave in the shock precursor, suggesting that the Fermi acceleration mechanism may switch off upon 
interaction with a molecular cloud. 

5. Discussion 

It has been s uggested th at the non-resonant cosmic ray current driven instability may play a significant role in the acceleration of 
cosmic rays (lBelill2004l) . In the case of strongly driven unstable modes, the mean field is unable to quench the growth of the waves 
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at 5B < Bq. The maximum achievable energy may be pushed beyond the Lagage-Cesarsky limit dLagage & Cesarskvl Il983h . not 
only by producing nonlinear turbulence to scatter the particles but also through amplification of the ambient field. 

In this paper we have shown that the non-resonant instability can have a very rapid growth rate under most young SNR condi- 
tions, in the interstellar medium. Thermal eff'ects are of little importance for nonrelativistic shocks, and large neutral fractions are 
required to prevent growth beyond the linear regime. However, efficient injection and acceleration are essential for this mechanism 
to occur. The largest uncertainty in our models arises in determining the strength of the driving term ^, which in our models depends 
only on ionisation fraction and not explicitly on the ion density. As a result the instability increases in more dense regions. It should 
be noted however, that energy losses close to the injection momentum can reduce cosmic ray production in over dense regions 
dPrurvet al.Lll996l) . 

The nonlinear evolution of field amplification via this instability requires further study. Current MHD simulations have only 
been performed with constant cosmic ray currents (Bell, 2004, 20 051), using periodic b oundary conditions. The first results of PIC 
simulations of this mechanism have recently been published (Niemiec & Pohll |2007|) . The dynamical range over which the non- 
resonant instability operates is dependent on the mass differences between magnetised and unmagnetised particles. Realistic proton 
electron mass ratios are essential in simulating this process. This is most likely the reason for the relatively low saturation value of 
6B ~ Bq. As the field becomes highly nonlinear, 6B > Bq, the cosmic ray trajectories will be deflected. Future PIC simulations will 
give further insight the evolution and saturation of the no n-resonant instability. 

Both the numerical MHD simulations of BelT (2004, 2005) and the PIC simulations of Nie miec & Pohll (|2007|) show that in the 
initial stages of the instability's development cavities and filaments are created along the direction of the mean magnetic field. It is 
cl ear from equa t ion (flOl i that the resulting cavities would be mostly populated by neutral species. This is in contrast with the results 
of lTagper et 311 ( 119951) 7 where the wave damping expels the neutrals from ionised flux tubes, although the waves they considered 
were not d riven by any mech anism. In highly dense, low ionisation regions, where the instability cannot exceed the damping, the 
analysis of lTagger et al.l d 1 9951) applies. 

The connection betwee n magnetic field amplification and the broadening of the narrow component Ha line has been discussed 
bv lGhavamian et a n (I2007h . These authors suggest that an observed broad component, of width 6v, such that (Jv/va > 1, implies 
a turbulent field component 6B/B > 1. We have demonstrated that the cosmic ray current driven instability is still capable of 
amplifying fields in the Balmer dominated shocks. This may provide a possible observational test of the amplification of the turbulent 
magnetic component in supernova shocks. 
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